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Abstract With a sky brightness monitor in Xinglong station of National Astronomical 
Observatories of China (NAOC), we collected data from 22 dark clear nights and 90 lu- 
nar nights. We first measured the sky brightness variation with time in dark nights, found a 
clear correlation between the sky brightness and human activity. Then with a modified sky 
brightness model of moon night and data from moon night, we derived the typical value for 
several important parameters in the model. With these results, we calculated the sky bright- 
ness distribution under a given moon condition for Xinglong station. Furthermore, we sim- 
ulated the moon night sky brightness distribution in a 5 degree field of view telescope (such 
as LAMOST). These simulations will be helpful to determine the magnitude limit, exposure 
time as well as the survey design for LAMOST at lunar night. 
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1 INTRODUCTION 

The night sky brightness of the observatory, which is not completely dark, is the major factor that constrains 
the limiting magnitude of a telescope and the exposure time given the demanded signal to noise ratio of the 
targets. For a fiber spectral instrument, sky subtraction is the crucial step in fiber spectra data reduction, 
usually this step relies on how the fiber sample the background which is generally thought to be homoge- 
nous within 1 degree in dark nights. Once there is a gradient in sky background (such as in moon night), 
much more effort will be needed in both designing the sky sampling fiber and data reduction. Thus a good 
estimation of sky background distribution is very important in designing a fiber spectroscopic survey with 
large field of view such as Large sky Area Multi-Object fiber Spectroscopic Telescope (LAMOST a.k.a. 
GuoShouJing telescope, Cui et al. (2012)). 

There are several sources that contribute to the night sky light after astronomical twilight, namely air- 
glow, scattering of starlight and zodiacal light, natural and artificial light pollution. For an astronomical 
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observatory the artificial light pollution is generally very small and the most dominative light pollution 
source is the moonlight when the moon is above the horizontal. Many previous work has been done to study 
the night sky brightness, for example: Sanchez et al. (2007),Schneebergeret al. (1979),Neugent & Massey 
(2010). However, most of the papers focus on the dark night sky brightness, only a few of them study 
the sky brightness in lunar nights. For LAMOST site, the Xinglong Station of the National Astronomical 
Observatories in China (NAOC), which is located 100 km northwest to Beijing at a longitude of 7h50ml8s 
east, a latitude of 40°23'36" north, and an altitude of 950 m, Liu et al. (2003) and Yao et al. (2012) had 
studied the night sky brightness , seeing, extinction, and available observational hours as well as their sea- 
sonal changes using photometric data of the BATC telescope. For the moon night sky, they use simple 
empirical model considering the moon phase, height and angular distance between the moon and target sky 
to eliminate the effect of the moon. The model used in their paper is too simple to estimate the detail sky 
brightness distribution in a 5 degree field of view ( as LAMOST). Krisciunas & Schaefer (1991) derived 
a more complicated model considering both the Rayleigh and Mie scattering of the moon light, two more 
variables, zenith distance of the moon and zenith distance of the sky were considered in their model, the 
accuracy reported in their paper was 8-23%. In the present work, with a modified moon night sky brightness 
model of Krisciunas & Schaefer (1991) and data from the measurement of a Sky Quality Meter (SQM) in 
Xinglong station, we derived the typical value for several important parameters in the model. With these 
results, we calculated the sky brightness distribution under a given moon condition. Furthermore, we use 
these results to predict the typical brightness gradient for LAMOST telescope 5 degree field of view. 

In section 2 we describe the sky brightness model. In section 3 the process about how we derive the 
brightness measurement data from the SQM is explained. In section 4 we present the fitting of parameters, 
including the dark zenith sky brightness V zen , the extinction coefficient k, the scattering coefficient PA 
and PB. In section 5 we calculate the sky brightness distribution when the moon is at a specific position 
and phase angle, using the typical value of deduced sky parameter . Furthermore the brightness variation 
distribution in a 5 degree field of view caused by the moon is predicted and the affection with the LAMOST 
is discussed. 

2 CALCULATION MODEL 

According to Krisciunas & Schaefer (1991), The sky brightness in lunar night can be divided into two part, 
the dark time night sky: 


B 0 (B zen , k, z) = B zen w-o *H x -i) X (1) 

and the contribution from the scattered moon light: 

B moon = f(p)I* 10-°- 4fe ^ (Z -)(l - IQ-0 ^X(Z)^ ( 2 ) 

where B zen is the dark time sky brightness at zenith in nanoLamberts (nL) which can be converted to 
magnitude V zen with equation (27) in Garstang (1989) , k is the extinction coefficient, f (p) is the scattering 
function at scattering angle p, I* is the brightness of the moon outside the atmosphere, can be expressed as 
a function of moon phase angle a : 
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Z is the zenith distance, X is the airmass, X. m ( Z m ) and X(Z) denote the airmass for the moon and the target 
sky respectively. The airmass X can be expressed as in Garstang (1989) : 

X(Z) = (1 - 0.96 sin 2 Z)~ 0 - 5 (4) 

The scattering function is composed of two types of scattering in the atmosphere, the Rayleigh scattering 
from atmospheric gases and the Mie scattering by atmospheric aerosols. The Mie scattering may degenerate 
into Rayleigh scattering when the size of atmospheric aerosols decreases. Instead, when the size increase 
the scattering turn out to be geometrical optics. The scattering function is proportional to the fraction of 
incident light scattered into a unit solid angle with a scattering angle, and also varies with different wave- 
length. According to Krisciunas & Schaefer (1991) and Chakraborty et al. (2005)), we define the scattering 
functions as 


f R (p) = 10 5 ' 36 (1.06 + cos 2 p) 

(5) 

/m{p) = 10 615_ «,wj hen(p > 10) 

(6) 

/m(p) = 6.2 x 10 7 p~ 2 , when(p < 10) 

(7) 

f(p) = PAxf M (p) + PBxf R (p) 

(8) 


In these functions, f R (p) and /m(p) are the Rayleigh and Mie scattering functions respectively, p is the 
scattering angle defined as the angular separation between the moon and the sky position. In equation 8, PA 
is the Mie scattering scale factor and PB is the Rayleigh scattering scale factor, they are proportional to the 
density of the Mie and Rayleigh scattering particles in the atmosphere. In Krisciunas & Schaefer (1991), 
the scattering function(equation 16 in that paper) have absorbed constant factor relating to unit conversions 
and normalizations. Since the particle densities may change with time and weather conditions from site to 
site, we use PA and PB as scale factors to scale the relative particle density to the site (Mauna Kea) in that 
paper, and they were treated as free parameters rather than fixed values in the data fitting program in section 
4 to deal with different weather conditions. If the scattering angle is small, the moonlight may directly enter 
the measurement instrument and equation 5 is not applicative in that case. 

Finally, the sky brightness in the moon night can be simply expressed as 

B = Bq + Bm 00n (9) 

Here, Bq is the dark night sky brightness in equation I and B moon is the brightness caused by the moonlight 
in equation 2. Combining all the equations before, the sky brightness can be expressed as 

B B(B zen , k , Z moon , Z sky , p, a, PA, PB) (10) 

Here B zen is the dark zenith sky brightness in nL; k is the extinction coefficient; Z moon is the moon zenith 
distance; Z sky is the zenith distance of sky position; p is the scattering angle; a is the moon phase angle; 
PA is the Mie scattering scale factor; PB is the Rayleigh scattering scale factor. 

The moon longitude and latitude can be calculated using Chapront ELP-2000/82, (Meeus (1991)). The 
equator coordinates (RA, DEC) can then be obtained. For application in this paper, the accuracy of the mean 
position of the moon is enough. The Greenwich Mean Sidereal Time can be calculated as 
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Here T is the Julian centuries from J2000.0. Then the hour angle of the moon can be calculated as 

Hi OC ai = CM ST — RA — LON (12) 

Here the geographic longitude (LON) is 117.575703° and the geographic latitude (LAT) is 40.3933333°. 
The Azimuth angle (A) and the elevation angle (h) of the moon can then be calculated as 

sin h = sin LAT sin DEC + cos LAT cos DEC cos Hi oca i 
cos hcosA = cos LAT sin DEC — sin LAT cos DEC cos Hi oca i (13) 

cos hsinA = — cos DEC sin Hi oca i 

3 SKY BRIGHTNESS DATA IN XINGLONG 

Sky Quality Meter (SQM, see http://www.unihedron.com/projects/darksky/) is a handy tool developed by 
Unihedron company to measure the sky brightness of visual light in mag/arcsec 2 . There is a near-infrared 
blocking fiber to contain the light to visual band, the transmission curve of SQM can be found in Fig22 
of Cinzano (2005)(http://www.lightpollution.it/download/sqmreport.pdf), the wavelength response is very 
broad, with half maximum of sensitive curve from 4000 to 6000A and peaked at 5400A. Due to the broad- 
ness of the response, the conversion to Johnson V band will depend on the spectral type. But for dark 
night or moon light sky brightness, the SQM magnitude is similar to V band, with an error of about 0. 1 
mag(Cinzano (2005)). Each SQM is calibrated by a NIST-traceable light meter, the absolute precision of 
each meter is believed to be 10% ( OAmag / ascsec 2 ). The difference in zero point between each calibrated 
SQM is also 10% (OAmag /ascsec 2 ). Add the above errors together, the overall sky brightness measured 
in V band is about 0.2 mag/arcsec 2 . 

A SQM is installed in the Xinglong station of NAOC. It measured the sky brightness in magnitude 
per square arcsecond per 6 minutes for more than a year. The Full Width Half Maximum (FWHM) of the 
angular sensitivity is 20°. The SQM is fixed in a metal framework and will not change its direction until 
reinstalled manually. Table 1 shows the direction of SQM in degrees during the days when the data was 
produced. 

Table 1 Direction of the SQM 


Start Date End Date 

Azimuth(°) Zenith Distance(°) 

2010-12-16 2011-12-14 1 80(South ) 

15 

2011-12-15 2012-01-12 

180(South) 

30 

2012-01-13 2012-02-13 

90(East) 

30 

2012-02-14 2012-03-12 

O(North) 

30 


Notes: The azimuth zero point is the north and the clockwise direction is positive. 


To measure the sky brightness of dark nights, we pick out data in the nights around new moon when the 
moon is below horizon, then we look up the observation log of LAMOST to reject those nights that was 
not marked as clear night, this give us about 22 nights. It will be interesting to explain how those data was 
influenced by light pollution from nearby cities and solar activities, but as we only have a handful nights for 
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Record time/hour 


Fig. 1 The average of sky brightness with night time from 22 possible clear dark nights. The 
vertical box at each data point shows 1 a error bar, vertical solid lines indicate the maximum 
and minimum value of that data point. There is a clear tendency that the night sky brightness is 
brighter in the first half of night, it goes darker in the second half due to less human activity . 


detail pattern of dark sky brightness. We still need time to accumulate more data to find the nightly pattern 
with seasons and directions. So, here we only show the average dark night brightness with time. We first 
correct the sky brightness to the zenith with equation I , then we set midnight as fiducial time, the data from 
different nights were averaged every 10 minutes from the fiducial. The results is plotted in figure I, there 
is a clear tendency that the sky brightness is brighter in the first half of the night, after midnight the sky 
brightness gets darker about 0.3 magnitude. This may be explained as the human activities decrease after 
midnight, then the light pollution is reduced at second half of night. The overall dark night sky brightness is 
about 21.6 ± 0.2 mag / arcsec 2 . The error bar for each data point in figure ! should reflect both the stability 
of the instrument and the nightly change of the sky condition, comparing with the sky brightness scatter in 
order of 0.4 magnitude in Liu et al. (2003) paper, also considering that our data were obtained in a period 
of about one year, the dominating contributor to this error bar should be the sky itself, the stability of SQM 
should be much smaller than 0.2 mag. 

The direction of the SQM was adjusted with simple tools so the accuracy is in several degrees. Those 
data were taken in more than 400 days with all kinds of weather conditions such as heavy cloudy days also 
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scattered moonlight, we need to exclude those situation which the scattering model can not work properly , 
we use several rules to filter the proper data set: 

1. Exclude the daytime records according to the astronomical twilight. 

2. Exclude the records when the altitude of the moon is less than 5°. 

3. Exclude the records when the angular distance is smaller than 32.5°, or the sky brightness is brighter 
than 18 {mag / arcs ec 2 ) to avoid direct incidence of moon light. 

4. Exclude the observation night whose record count is smaller than 25 to avoid night with heavy cloud. 

5. Exclude the observation night whose observed brightness curves show irregular transition or large dis- 
persion than 30% to avoid drastic weather change or partial cloud blocking. 

In the end, we obtained sky brightness data for 114 moon nights from Dec 16th 2010 to Mar 12th 2012. 
Each observation night contains about 42 records. 


4 PARAMETERS FITTING 


This following work is mainly based on equation 10. Here the zenith distance of the moon Z moon , the 
zenith distance of the sky position Z s k y , the angular separation between the moon and the sky position p, 
and the moon phase angle a in degree are managed as input arguments, since they can be calculated from the 
observatory geographical location and time. The dark zenith sky brightness V zen , the extinction coefficient 
k, the Mie scattering scale factor PA and the Rayleigh scattering scale factor PB are treated as parameters 
to be determined. The output result is sky brightness in nanoLamberts, which is provided through the SQM 
measurements. 

In order to solve the non-linear least squares problem, we collect the observed sky brightness records 
on each observation night as one data set. The number of data records (N) in each data set is approximately 
42. Using the observatory geographical location and time information, the moon and sky position related 
input values can be estimated and attached to each record. It is assumed that for one night the parameters 
( V zen , k, PA and PB) do not change very much so they are treated as constants. We can determine the value 
of these parameters for each night that will give a minimum value of the squared 2-norm residual, which is 
defined as 


i=N 

RESNORM = {ModelBi - ObservedBi ) 2 (14) 

t=i 

The initial value of the parameters ( V zen , k, PA and PB) are set as 21.4, 0.23, 1 and 1 respectively. The 
zenith magnitude is limited between 0 and 25; the extinction coefficient is limited between 0.01 and 8; the 
parameter PA and PB are limited between 0 and 25. Finally we obtained the fitting results for 114 data sets. 
If we reject data sets with returned parameters at or close to the limitation boundary, 90 data sets remained. 

Figure 2 shows the model brightness against the measured brightness for the data records of all the data 
sets from 1 14 nights, containing more than 5000 records. We define the Relative Fitting Variation (REV) as 


RFV 


\ 


7— AT . „ 

1 / ModelBi — ObservedBi 

N l ObservedBi 

2=1 X 


2 


(15) 


The total relative fitting variation defined by equation (15) in Figure 2 is 12%. The RFVs for each data set 


Moon night sky brightness simulation for Xinglong station 


7 



Fig. 2 The model brightness against the observed brightness. About 5200 data records from the 
1 14 nights are presented. The total relative fitting variation defined by equation (15) is 12%. 


results have a RFV smaller than 5%. That means for most of the remaining nights, the model brightness 
matches the observed brightness with a nice accuracy. 


c 
3 

o 
(_> 



Relative Fitting Variation 


Fig. 3 Histogram of relative fitting variation of each observation night. Only the the last remain- 
ing 90 nights are presented. Most of the variations are smaller than 5%. 


Figure 4 shows examples of fitting curve from four nights. The moon phase angle here is expressed 
as the Sun-Moon-Earth angle in degrees. The relative fitting variations are 2.9%, 4.7%, 3.5% and 10% 
respectively, as marked in the figure. It is obvious that the bottom-right figure has the maximum brightness 
because the phase angle is the smallest among the four figures, while the top-right one has the minimum 
brightness. The major influence factor of the curve trend is the moon/sky angular separation. As we can see 
from these figures that when the moon moves close to the sky position, the sky background becomes bright. 
Since the sky direction and the moon declination for one night is basically fixed, the angular separation is 
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of view of a quasi-meridian telescope, such as LAMOST, will be very bright. The results shows that this 
fitting model works well for various phase angle and location of the moon. 

Nights with bad fitting results are due to following reasons: 

1 . The weather condition. As this measurement instrument keeps working without interruption for years, 
the data is produced on days under various kinds of weather condition. These affection cannot be cor- 
rected using the presented models. 

2. The large angular sensitivity of the SQM. Within that large field of view the light from stars and the 
planets, especially the moonlight when the moon is close, will greatly affect the measure results. As a 
result, the measured sky brightness shows a sharp increase as the moon moves close to the sky position. 

3. The reflection of the moonlight from the surrounding buildings. 

4. The light pollution from the nearby cities. The sky brightness of Xinglong station is mainly polluted 
by the city lights from Beijing, Xinglong, and Chengde. As we can see from figure I , the artificial light 
pollution may change the sky brightness about 0.3 mag at night. Those affections are ignored in the 
model since the sky brightness we are considering is at least one magnitude brighter than the dark night. 
The sky brightness in the absence of moonlight is expressed as equation (1). 

After iterative fitting, the result histogram of the estimated parameters including the dark zenith sky 
brightness V zen , extinction coefficient k, scattering coefficients PA and PB are shown in Figure 5. Only the 
last remaining 90 data sets are presented. The dark zenith sky brightness, V zen , presents a typical value 
of 21.4 mag / arcsec 2 . This result agrees with the zenith sky brightness (21.6 ± 0.2mag / arcsec 2 ) in the 
winter we measured in section 3 from dark nights. This also basically fit the sky brightness in the V band 
obtained from BATC Polaris monitoring data (Yao et al. (2012)), which is about 21 mag / arcsec 2 at the 
North Pole. Correcting the airmass using equation 1 for 0.36 magnitude, the results agree very well. The 
extinction coefficient, k, presents a typical value of 0.23. This is in the acceptable range according to the 
measurement results in Xinglong (Liu et al. (2003); Yao et al. (2012)). The Mie scattering scale factor, PA, 
presents a typical value of 1.50. The Rayleigh scattering scale factor, PB, presents a typical value of 0.90. 

Figure 5 shows the estimated parameters as a function of date from January 1 2011. The data were 
collected within one year, there is an obvious lack of data in the summer due the bad weather. Since there 
are not enough data, it’s hard to tell the system change of those parameters with seasons. However we can 
still find some indication that the dark zenith sky is darker in the winter than in the summer and the sky 
transparency in the autumn and winter is better than those in the spring and summer, which is also consistent 
with previous work, e.g. Liu et al. (2003) and Yao et al. (2012). While the scattering scale factor PA and PB 
show no obvious seasonal variation. The factor PA has a larger dispersion than PB. This is mainly because 
the density of the Mie scattering particles (i.e. dust) has a bigger variation than the Rayleigh molecule. More 
data are need to better resolve the seasonal change of these parameters. 

In Liu et al. (2003) paper, they found a linear relation between the sky brightness and extinction coef- 
ficient k. While in our results( see figure 7), the correlation is not found. Since the data in Liu et al. (2003) 
paper were collected around Polaris, there is a fix zenith distance about 50° in Xinglong station, consid- 
ering equation 1, we conclude that the linear correlation agrees with the equation while the measured sky 
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Model Brightness 

— □— Observed Brightness 



Fig. 4 Sky brightness fitting curves. The X axis is the serial number for the measured data point. 
Related parameters, including observation date, moon age, moon phase angle in degrees, moon 
declination and relative fitting variation RFV, are printed in the figures. The definition of RFV is 
described in equation ( 15). The SQM directions are listed in table 1 . 
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Moonless Zenith Sky Brightness (magnitude) 



Extinction coefficient, k 



02468 10 12 14 16 0 1 234567 


The Mie Scattering Ratio Factor, PA The Rayleigh Scattering Ratio Factor, PB 


Fig. 5 Histogram of the estimated parameters including the dark zenith sky brightness V zen , 
extinction coefficient k, scattering coefficients PA and PB. Only the last remaining 90 data sets 
are presented. The typical value for these parameters can be read from the figures as: V zen =2l.4, 
k=0.23, PA=1.5 and PB=0.90 
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Days from UTC 201 101 01 12:00:00 Days from UTC 20110101 12:00:00 




Days from UTC 20110101 12:00:00 Days from UTC 20110101 12:00:00 


Fig. 6 Estimated parameters on days. Top left: zenith sky brightness vs. time. Top right: extinc- 
tion vs. time. Bottom left: Mie scattering factor PA vs. time. Bottom right: Rayleigh scattering 
factor PB vs. time. Only the last remaining 90 data sets are presented. The days are calculated 
from 1 January 2011, 12:00:00. Days close to 0 or 365.25 are in the winter. Days near 180 are in 
the summer. 



Fig. 7 extinction coefficient k vs. zenith sky brightness 14^. 


12 


S. Yao, et. al. 


5 SKY BRIGHTNESS DISTRIBUTION 

5.1 Sky brightness estimation 

Estimation of sky background is important for a telescope working at lunar night. The background bright- 
ness is one of the key factor to determine the limiting magnitude and exposure time. For a survey telescope 
working at lunar night, knowing the background distribution will help astronomers to design how faraway 
from the moon should the telescope point to in order to reduce the influence of the moon. By applying 
the typical parameters of fitting results to equation (10), we can calculate the sky brightness for Xinglong 
station given sky position and time within lunar night, except the region where the moon is very close to 
the sky position (< 10°). Here we give two examples in table 2 and table 3. In table 2 the moon is assumed 
to be with a phase angle of 20°, an hour angle of -30° and a declination of -10°. The brightness distribution 
in a large region with hour angle from -35° to 35° and declination from -10° to 70° is presented. For an in- 
tuitive view, sky brightness distribution for moon phase angle 20°, declination —20° and hour angle —30° 
is plotted in figure 8. Table 3 is for a moon on the meridian with phase angle 90°, declination 20°. The 
position within 10° of the moon was left blank in the table. The brightness in table 2 and 3 will be used to 
set the current magnitude limit and estimate exposure time for FAMOST survey in each lunar night. They 
could also help to determine where the telescope should point to given the scientific request. For example, 
to reach magnitude 17 in bright night, the background should be no brighter than 19 mag / arcsec 2 . In table 
2 and figure 8, the moon is close to full moon, sky should be at least 65° away to reach 19 mag j arcsec 2 , 
comparing with table 3 for which the moon is at half moon, the sky is as deep as 19 mag / arcsec 2 when 
15° from the moon. From table 2, 3 and figure 8, we can see the brightness gradient is larger when closer to 
the moon. The sky brightness drops faster in declination directions than hour angle direction. In declination 
direction, sky brightness drops faster when pointing to the zenith than to the horizon. 



Fig. 8 Distribution of sky brightness in magnitude, the input parameters are: V zen =2 1 .4, K=0.23, 
Moon phase=20°, DEC=-20°, hour angle=-30°, PA=1.5, PB=0.90 . 
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Table 2 The sky brightness distribution in Xinglong station when the moon is 2 
hours before the transit. The brightness is in V band mag/arcsec 2 . input conditions 
are listed below the table. 


Dec\h 

-35 

-30 

-25 

-20 

-15 

-10 

-5 

0 

5 

10 

15 

20 

25 

30 

35 

-10 


Moon 



17.16 

17.39 

17.59 

17.77 

17.93 

18.08 

18.20 

18.30 

18.38 

18.45 

18.49 

-5 




17.06 

17.29 

17.50 

17.70 

17.88 

18.03 

18.17 

18.29 

18.40 

18.48 

18.55 

18.60 

0 

17.02 


17.10 

17.27 

17.45 

17.64 

17.83 

17.99 

18.14 

18.28 

18.39 

18.49 

18.58 

18.64 

18.69 

5 

17.29 

17.30 

17.37 

17.50 

17.65 

17.81 

17.97 

18.12 

18.26 

18.38 

18.49 

18.59 

18.67 

18.73 

18.78 

10 

17.55 

17.57 

17.63 

17.72 

17.84 

17.98 

18.11 

18.25 

18.37 

18.49 

18.59 

18.68 

18.75 

18.81 

18.86 

15 

17.79 

17.81 

17.86 

17.94 

18.04 

18.15 

18.26 

18.38 

18.49 

18.59 

18.68 

18.76 

18.83 

18.89 

18.93 

20 

18.01 

18.03 

18.07 

18.14 

18.22 

18.31 

18.41 

18.51 

18.60 

18.69 

18.77 

18.85 

18.91 

18.96 

18.99 

25 

18.20 

18.23 

18.27 

18.32 

18.39 

18.46 

18.54 

18.63 

18.71 

18.79 

18.86 

18.92 

18.98 

19.02 

19.05 

30 

18.38 

18.40 

18.44 

18.48 

18.54 

18.60 

18.67 

18.74 

18.81 

18.88 

18.94 

18.99 

19.04 

19.08 

19.10 

35 

18.54 

18.56 

18.59 

18.63 

18.68 

18.73 

18.79 

18.85 

18.90 

18.96 

19.01 

19.06 

19.10 

19.13 

19.14 

40 

18.68 

18.70 

18.73 

18.76 

18.80 

18.84 

18.89 

18.94 

18.99 

19.03 

19.07 

19.11 

19.14 

19.17 

19.18 

45 

18.80 

18.82 

18.84 

18.87 

18.91 

18.94 

18.98 

19.02 

19.06 

19.10 

19.13 

19.16 

19.18 

19.20 

19.20 

50 

18.90 

18.92 

18.94 

18.97 

19.00 

19.03 

19.06 

19.09 

19.12 

19.15 

19.17 

19.19 

19.21 

19.22 

19.22 

55 

18.99 

19.01 

19.03 

19.05 

19.07 

19.09 

19.12 

19.14 

19.16 

19.18 

19.20 

19.22 

19.23 

19.23 

19.23 

60 

19.05 

19.07 

19.09 

19.11 

19.12 

19.14 

19.16 

19.18 

19.19 

19.21 

19.22 

19.23 

19.23 

19.23 

19.22 

65 

19.10 

19.11 

19.13 

19.14 

19.16 

19.17 

19.18 

19.20 

19.21 

19.21 

19.22 

19.22 

19.22 

19.22 

19.21 

70 

19.12 

19.14 

19.15 

19.16 

19.17 

19.18 

19.19 

19.20 

19.20 

19.20 

19.21 

19.21 

19.20 

19.20 

19.19 


Notes: Unit:mag. V zen =21.4, k=0.23, moon phase=20°, hour angle=-30°, DEC=-10°, PA=1.5, PB=0.90 


5.2 Brightness difference within a 5 degree field of view 

The Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) is an optical survey telescope 
with a field of view of 5 degrees. There are 16 spectrographs each holding 250 fibers, each spectrograph 
occupy a certain field of view( about 1°) on the LAMOST focal plane. In such a large field of view, the 
background gradient can not be ignored. In most fiber spectra observations, the sky is sampled by dedicated 
fibers, then the sky in the object fiber is subtracted using a spectrum composited from those sky fibers. This 
step works well when the sky is homogenous within the field of view of spectrograph, in the case of moon 
night, the sky background gradient should be considered in the data reduction step as well as in the step 
that designing how the sky fiber sampled the background. In this work, we try to estimate the brightness 
difference caused by the moonlight within the field of view of LAMOST, so as to provide reference basis for 
observation strategy decision as well as data reduction. From table 2 and 3, it is easy to tell that the larger the 
angular separation, the less the sky brightness gradient is. To calculate the maximum brightness difference 
within the field of view, we need to find out the point with the maximum and minimum brightness. The 
points with extreme value was found by a step by step search along the edge of field of view. We show an 
example of results in table 4. The input parameters are the same as in table 3 except that the moon phase 
angle is 20° rather than 90°. Each item in table 4 is the max difference in LAMOST 5 degree field of 
view. The max difference in table 4 is 0.36 mag/ 5°, so in the field of view of one LAMOST spectrograph 
(about 1°) , the gradient will be about 0.07 mag. Which means the gradient even in one spectrograph is 7%, 
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Table 3 The sky brightness distribution when moon phase angle is 90 degrees 


Dec\h 

-35 

-30 

-25 

-20 

-15 

-10 

-5 

0 

5 

10 

15 

20 

25 

30 

35 

-10 

19.53 

19.50 

19.46 

19.42 

19.38 

19.35 

19.33 

19.32 

19.33 

19.35 

19.38 

19.42 

19.46 

19.50 

19.53 

-5 

19.55 

19.50 

19.45 

19.39 

19.34 

19.29 

19.26 

19.24 

19.26 

19.29 

19.34 

19.39 

19.45 

19.50 

19.55 

0 

19.56 

19.50 

19.43 

19.35 

19.28 

19.21 

19.16 

19.14 

19.16 

19.21 

19.28 

19.35 

19.43 

19.50 

19.56 

5 

19.56 

19.49 

19.40 

19.31 

19.21 

19.11 

19.04 

19.01 

19.04 

19.11 

19.21 

19.31 

19.40 

19.49 

19.56 

10 

19.57 

19.48 

19.38 

19.26 

19.14 

19.01 

18.90 


18.90 

19.01 

19.14 

19.26 

19.38 

19.48 

19.57 

15 

19.58 

19.48 

19.37 

19.23 

19.08 

18.92 




18.92 

19.08 

19.23 

19.37 

19.48 

19.58 

20 

19.59 

19.49 

19.38 

19.24 

19.08 



Moon 



19.08 

19.24 

19.38 

19.49 

19.59 

25 

19.62 

19.52 

19.41 

19.28 

19.13 

18.97 




18.97 

19.13 

19.28 

19.41 

19.52 

19.62 

30 

19.66 

19.57 

19.47 

19.35 

19.23 

19.11 

19.02 


19.02 

19.11 

19.23 

19.35 

19.47 

19.57 

19.66 

35 

19.71 

19.63 

19.54 

19.45 

19.36 

19.27 

19.21 

19.19 

19.21 

19.27 

19.36 

19.45 

19.54 

19.63 

19.71 

40 

19.76 

19.70 

19.63 

19.56 

19.49 

19.43 

19.39 

19.38 

19.39 

19.43 

19.49 

19.56 

19.63 

19.70 

19.76 

45 

19.82 

19.77 

19.72 

19.67 

19.62 

19.58 

19.55 

19.55 

19.55 

19.58 

19.62 

19.67 

19.72 

19.77 

19.82 

50 

19.88 

19.84 

19.81 

19.77 

19.74 

19.71 

19.69 

19.69 

19.69 

19.71 

19.74 

19.77 

19.81 

19.84 

19.88 

55 

19.94 

19.91 

19.89 

19.86 

19.84 

19.82 

19.81 

19.81 

19.81 

19.82 

19.84 

19.86 

19.89 

19.91 

19.94 

60 

19.99 

19.97 

19.95 

19.94 

19.92 

19.91 

19.91 

19.90 

19.91 

19.91 

19.92 

19.94 

19.95 

19.97 

19.99 

65 

20.03 

20.02 

20.01 

20.00 

19.99 

19.99 

19.98 

19.98 

19.98 

19.99 

19.99 

20.00 

20.01 

20.02 

20.03 

70 

20.07 

20.06 

20.06 

20.05 

20.05 

20.04 

20.04 

20.04 

20.04 

20.04 

20.05 

20.05 

20.06 

20.06 

20.07 


Notes: Unit:mag. V ze „=2l A, k=0.23, moon phase=90°, hour angle=0°, DEC=20°, PA=1.5, PB=0.90 


difference in the table is about 0.05 mag which means about 1% gradient in one spectrograph. Generally, the 
sky subtraction accuracy in fiber spectra data reduction is larger than 2%, 1% difference can be acceptable in 
most cases. From the above discussion, it’s better to choose the position with smaller background gradient 
to alleviate data reduction difficulty in designing a survey like LAMOST. Figure 9 shows an example of sky 
brightness distribution inside LAMOST 5 degree field of view, we calculate the brightness for each fiber. 
The moon phase angle is 20°, the distance between the moon and the field center is about 30° as marked 
in the figure. The moon is 67° to the north east, as show by the lower right icon. As we can see, there is 
a brightness gradient about 0.16 magnitude, but the gradient direction is not the direction from the moon 
to the sky position. As we pointed out in section 5.1, this is because when the airmass increase against the 
zenith direction( i.e. north in this figure), the moon light gets more scattered by the atmosphere, this will 
change the direction of gradient a bit to the south . 


6 SUMMARIES 

We use a SQM to study the night sky brightness in Xinglong station. From the collected data from Dec 2010 
to Mar 2012, we selected 22 dark clear nights to study the sky brightness variation with time, we found a 
clear correlation between the dark night sky brightness with human activity. We also study the lunar sky 
brightness model of Krisciunas & Schaefer (1991), by modifying the relative scale factors of Rayleigh and 
Mie scattering respectively, we successfully fitted the sky brightness data of SQM in Xinglong station with a 
relative fitting variation of 12%. We estimate the related parameters in 90 nights. According to the results we 
can see that in this observatory the typical V band dark zenith sky brightness is about 21.4 mag/ascsec 2 \ 


Moon night sky brightness simulation for Xinglong station 
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Table 4 Brightness difference within LAMOST field of view 
when the moon is over the zenith 


Dec\h 

-35 

-30 

-25 

-20 

-15 

-10 

-5 

0 

5 

10 

15 

20 

25 

30 

35 

-10 

0.05 

0.06 

0.06 

0.07 

0.08 

0.09 

0.09 

0.09 

0.09 

0.09 

0.08 

0.07 

0.06 

0.06 

0.05 

-5 

0.06 

0.07 

0.08 

0.10 

0.10 

0.11 

0.11 

0.11 

0.11 

0.11 

0.10 

0.10 

0.08 

0.07 

0.06 

0 

0.08 

0.09 

0.10 

0.12 

0.13 

0.13 

0.14 

0.14 

0.14 

0.13 

0.13 

0.12 

0.10 

0.09 

0.08 

5 

0.09 

0.11 

0.12 

0.14 

0.15 

0.16 

0.16 

0.16 

0.16 

0.16 

0.15 

0.14 

0.12 

0.11 

0.09 

10 

0.11 

0.12 

0.14 

0.16 

0.17 

0.18 

0.23 


0.23 

0.18 

0.17 

0.16 

0.14 

0.12 

0.11 

15 

0.12 

0.14 

0.15 

0.17 

0.19 

0.29 




0.29 

0.19 

0.17 

0.15 

0.14 

0.12 

20 

0.13 

0.15 

0.17 

0.19 

0.20 



Moon 



0.20 

0.19 

0.17 

0.15 

0.13 

25 

0.13 

0.15 

0.17 

0.19 

0.21 

0.36 




0.36 

0.21 

0.19 

0.17 

0.15 

0.13 

30 

0.13 

0.15 

0.17 

0.19 

0.21 

0.23 

0.26 


0.26 

0.23 

0.21 

0.19 

0.17 

0.15 

0.13 

35 

0.13 

0.15 

0.17 

0.18 

0.20 

0.21 

0.22 

0.23 

0.22 

0.21 

0.20 

0.18 

0.17 

0.15 

0.13 

40 

0.13 

0.14 

0.16 

0.17 

0.18 

0.19 

0.20 

0.20 

0.20 

0.19 

0.18 

0.17 

0.16 

0.14 

0.13 

45 

0.12 

0.14 

0.15 

0.16 

0.17 

0.18 

0.18 

0.18 

0.18 

0.18 

0.17 

0.16 

0.15 

0.14 

0.12 

50 

0.11 

0.12 

0.13 

0.14 

0.15 

0.16 

0.16 

0.16 

0.16 

0.16 

0.15 

0.14 

0.13 

0.12 

0.11 

55 

0.10 

0.11 

0.12 

0.13 

0.13 

0.14 

0.14 

0.14 

0.14 

0.14 

0.13 

0.13 

0.12 

0.11 

0.10 

60 

0.09 

0.10 

0.11 

0.11 

0.11 

0.12 

0.12 

0.12 

0.12 

0.12 

0.11 

0.11 

0.11 

0.10 

0.09 

65 

0.08 

0.08 

0.09 

0.09 

0.10 

0.10 

0.10 

0.10 

0.10 

0.10 

0.10 

0.09 

0.09 

0.08 

0.08 

70 

0.07 

0.07 

0.07 

0.08 

0.08 

0.08 

0.08 

0.08 

0.08 

0.08 

0.08 

0.08 

0.07 

0.07 

0.07 


Notes: Unit:mag. V zm =l\ .4, k=0.23, moon phase=20°, hour angle=0°, DEC=20°, PA=1.5, PB=0.90 

scale factor is about 0.90. With the model and those typical parameters for Xinglong station, we could 
estimate the sky brightness distribution for any given time in moon night for the LAMOST site. We then 
calculated the gradient within the LAMOST 5 degree field of view. The result shows that sky brightness 
increases quickly as the distance of the moon is smaller. The increasing zenith distance will also enhance 
the brightness within a quantity much smaller than the influence of the lunar separation angle. These results 
will help in designing the LAMOST bright night survey, determine the location of sky fiber in the focal 
plane as well as data reduction in LAMOST survey. 
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Fig. 9 Simulation of brightness distribution in LAMOST focal plate. The moon is 2 hours before 
the zenith and the declination is 20°. The telescope is pointing to transit with declination of 10°. 
PA = 1.5, PB = 0.90, V zen = 21.4, k = 0.23, moon phase angle a= 20°. The fiber location is 
assumed to be placed in the sky according to their light path to make the figure more intuitive. 
Affected by the factor relating the zenith distance of the sky position, the moon direction does 
not match with the gradient direction of brightness. The increasing direction of brightness shifts 
from the moon to the horizon slightly. The variation becomes larger at the horizon. 
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